B-mode CMB Polarization from Patchy Screening during Reionization 
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_B-modes in CMB polarization from patchy reionization arise from two effects: generation of 
polarization from scattering of quadrupole moments by reionization bubbles, and fluctuations in the 
screening of i5-modes from recombination. The scattering contribution has been studied previously, 
but the screening contribution has not yet been calculated. We show that on scales smaller than the 
acoustic scale {£ > 300), the B-mode power from screening is larger than the B-mode power from 
scattering. The ratio approaches a constant ~2.5 below the damping scale {£ > 2000). On degree 
scales relevant for gravitational waves {£ < 100), screening _B-modes have a white noise tail and are 
subdominant to the scattering effect. These results are robust to uncertainties in the modeling of 
patchy reionization. 



I. INTRODUCTION 

In linear theory, the presence of a curl or i?-niode pat- 
tern in the polarization of the CMB is an indication that 
there are contributions from gravitational waves or vec- 
tor perturbations that, unlike density perturbations, can 
impart a sense of handedness to the polarization. Beyond 
linear theory, it is well known that second order effects 
such as gravitational lensing can produce i?-modes from 
density fluctuations [1]. In fact, any process that mod- 
ulates the polarization amplitude, direction, or position 
on the sky can generate B-modes from the intrinsic E- 
modes at recombination (e.g. [2]). 

In this Brief Report, we study one such contribution: 
the amplitude modulation due to patchy screening of the 
primary polarization due to inhomogencous reionization. 
Whereas the analogous effect of patchy generation of po- 
larization from scattering of the quadrupole anisotropy 
during reionization has been well-studied in the litera- 
ture [3-8] , patchy screening is usually neglected as a small 
contribution. While negligible at large angles, we show 
here that the two effects are always comparable in magni- 
tude beyond the damping tail with ~ 2.5 times as much 
power in the screening effect. 



II. PATCHY SCREENING 

A. Formalism 

Scattering of CMB radiation during reionization out of 
the line of sight suppresses the primary temperature and 
polarization anisotropy from recombination as where 
T is the Thomson optical depth. If r varies across the line 
of sight this suppression itself introduces anisotropy as 
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where T is the temperature fluctuation and Q and U are 
the polarization Stokes parameters, and we have assumed 
that the dominant temperature fluctuations from recom- 
bination are below the angular scale subtended by the 
horizon during reionization. These recombination fluctu- 
ations can be decomposed in multipole moments in the 
usual way with the additional assumption that the po- 
larization contains £'-modcs only 
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Decomposing the anisotropic part of the optical depth 
into multipole moments 
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we obtain in the (T(n) — r) ^ 1 limit 
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where eivv = [1 + (-l)^+^'+^"]/2 and ou-v = [1 
(— 1)^+^ +^ ]/2 picks out even and odd triplets and 



Here we have given the contribution to the anisotropy 
from the £ > terms in e""^: 
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where X G {T, B}. The additional power spectra con- 
tributions 
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We evaluate the recombination power spectra by setting 
T = 0. These sums can be efficiently evaluated in position 
space as shown in the Appendix. 

Note that in the decomposition defined in Eqn. (6), 
the isotropic screening term gains a contribution from 
the rms fluctuations in t 
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This term is ordinarily not included in the standard 
calculation but is a correction that is second order in 
(r(n)-f). 



B. B-Mode Scaling 

We now focus on the B-mode generation from screen- 
ing and its relationship to other secondary effects. To 
extract the large and small angle scaling behavior of the 
_B-modes it is useful to take the fiat-sky approximation 
where Fourier moments replace harmonic coefficients 
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where 0£ denotes the angle between i and the axis on 
which Q is defined. The power spectrum 
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then becomes for the screening modes 
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where we have set 4>i = for convenience. 

The i3-mode power spectrum therefore represents a 
convolution of recombination and r power spectra. Its 
main properties are defined by the well-determined shape 
of the primary .E-mode power spectrum. In particular, 
the primary i?-modes have little power above the acous- 
tic scale £ < £a ^ 300 and below the damping scale 
i > £d ^ 2000. In these two limits, the screening B- 
modes take on a particularly simple form. For £ ^ £a, 
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is white noise in form. Note that this property is in- 
dependent of the amount of low t power in CJ"^. Only 
optical depth fluctuations on scales comparable to £a can 
modulate the polarization into large angle i?-modes. 

For £ ^ £0, the B- modes reflect the shape of the t 
power spectrum 
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with an amplitude that is determined by the total power 
in the primary iJ-modes 
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Under the damping tail, the primary i?-modes are much 
smoother than r fluctuations and so i^rms simply repre- 
sents the typical level of i?-modc that is modulated into 
small scale _B-mode polarization by r. 

It is useful to compare these scalings with those of 
the two other well-known sources of secondary i?-modes: 
gravitational Icnsing and patchy scattering of quadrupole 
anisotropy at reionization. Lensing obeys a similar form 
with deflection angles playing the role of t and polariza- 
tion gradients playing the role of E. In particular, it has 
a white noise form for £ <^ £a and follows the deflection 
power spectrum for £ ^ £d with an amplitude set by the 
total gradient power in the primary polarization [1, 9]. 

For patchy scattering of the quadrupole fluctuations 
during reionization [3] 
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for scales much smaller than the reionization bump £ 3> 
£ii « 20. Here Qrms is the rms quadrupole at reionization 
and Teff is the typical optical depth between z = and 
the scatterers, so that Toff < f . The crucial difference 
between this effect and the patchy screening effect is that 
the quadrupole is coherent on large angles <C £a and 
so the i?-modcs are in the small scale limit for all relevant 
scales. 
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In particular for ^ ^ 
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For the WMAP5 cosmology [10], VLbK^ = 0.02265, 
^mh'^ = 0.137, h = 0.701, 17a 0.721, f = 0.084, 
As = 2.16 X 10^^, Us = 0.96, the rms fluctuations are 
-E'rms = 6.9/iK and Qrms = 17.9^K and so the ratio of 
power approaches a factor of ~ 2.5. Independently of 
the form of CJ^ , the screening B-modes dominate the 
scattering ones on scales near the damping tail and be- 
low. 

Conversely, independently of the power in CJ'^ at mul- 
tipoles £ < the screening B-modes fall as white noise 
whereas the scattering ones can continue to rise if there is 
large-scale power in r. Therefore only the scattering B- 
modes are relevant as contamination to the gravitational 
wave signal at I ^ 100. 

We show examples of these behaviors in Fig. 1. We 
take two representative models, a "fiducial" model in 
which CJ'^ has only small scale power contributed by rel- 
atively small ionization patches, and a "maximal" model 
which has substantial power at £ ~ 100 from patches that 
are substantially larger than those expected in current 
models of reionization. The large-bubble model maxi- 
mizes the scattering i?-modc contamination to the grav- 
itational wave signal [7] . 

More precisely, in both models the ionization history 
Xe{z) is the same as the fiducial model in [11], where the 
parameter Aj, controls the duration of partial ionization. 
Specifically, our choice of /S.y = 19 gives a range of 6 < 
z < 14. The total optical depth to recombination is taken 
to be f = 0.084. During partial reionization, the ionized 
regions are represented by spherical bubbles with a log- 
normal distribution [12, 13] with characteristic size given 
by R (in Mpc) and distribution width given by ax^R. We 
take {-R, (Tini?J = {5jlii(2)} in the fiducial model, and 
{-R, (Tinfl} = {30,ln(2)} in the large-bubble model. We 
assume that the bubbles are linearly biased tracers of the 
large-scale matter density field, with bias 6 = 6 using the 
construction of Ref. [14]. Note that rms fluctuations in 
T for the small and large-bubble models are Trms ~ 0.007 
and Trms = 0.017 and so both satisfy the condition that 

'^rms 

For reference, the small scale limits to the other power 
spectra for patchy screening can be obtained by replacing 
with Tj^g for temperature in Eqn. (14) where 
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whereas the E polarization and B polarization power are 
equal. The TE cross power spectrum and correlation 
coefficient is zero in this limit. 
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FIG. 1: B-mode power spectra. Three different contri- 
butions are being plotted; lensed B-modes, patchy B- 
modes from Thomson scattering, and patchy _B-modes from 
screening. The model is the WMAP5 cosmology with 
(a) the fiducial reionization model of {r, Ay, 6, _R, (Jin h} = 
{0.084, 19.0, 6.0, 5 Mpc, ln(2)}^and (b) a maximal model with 
the same parameters except -R = 30 Mpc. For reference the 
T power spectrum in dimensionless units is also plotted. 



III. DISCUSSION 

Whereas the patchiness of reionization itself re- 
mains theoretically uncertain and observationally uncon- 
strained, its effects on CMB polarization are well defined. 
In addition to the well-known patchy scattering of the 
quadrupole moment, the patchy screening of primary po- 
larization produces a i?-mode. 

Independently of the form of the patchiness, this con- 
tribution is larger than the scattering effect by a factor 
of ~ 2.5 in power on scales below the damping tail. Con- 
versely, on scales above the acoustic scale, the patchy 
screening effect always falls as white noise unlike the 
scattering effect. Consequently screening _B-modes are 
relatively more important when considering reionization 
contributions at arcminute scales but are largely unim- 
portant for the degree scales relevant for gravitational 
wave studies. 
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APPENDIX A: EFFICIENT POSITION-SPACE 
FORMS 

The harmonic-space forms for the power spectra given 
in Eqn. (8) naively require computational cost 0{£^^^) to 
compute the power spectrum for all multipoles £ < ^max- 
For the calculations in this paper (with €max ~ 10"*), 
we have found it convenient to use position-space forms 
which are mathematically equivalent but have reduced 
computational cost O(^max)- 
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where the correlation functions are defined by: 



An 



= ±^[^)cr(-^'4i,(0), (A2) 



where (i2±2i ^o±2 ^^'^ ^oo ^^'■^ reduced Wigner D- 
functions. The integrals in Eqn. (Al) can be done ex- 
actly, using Gauss-Legendre quadrature with [(Si'max + 
l)/2] points. 
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